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Abstract 

We construct static axially symmetric solutions of SU(2) Einstein- Yang-Mills- 
dilaton theory. Like their spherically symmetric counterparts, these solutions 
are nonsingular and asymptotically flat. The solutions are characterized by the 
winding number n and the node number k of the gauge field functions. For fixed 
n with increasing k the solutions tend to "extremal" Einstein-Maxwell-dilaton 
black holes with n units of magnetic charge. 
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1 Introduction 



SU(2) Einstein- Yang-Mills (EYM) theory possesses regular static spherically symmetric 
solutions fl| . These solutions are aymptotically flat and have non-trivial magnetic gauge 
field configurations, but no global charge. They have a high-density interior region, 
followed by a near-field region with approximately Reissner-Nordstr0m metric and a 
far- field region with approximately Schwarzschild metric [IJ. 

To every regular solution in SU(2) EYM theory, there exists a corresponding family 
of black hole solutions with regular event horizon xh > [0. Outside their event 
horizon these black hole solutions possess non-abelian hair. Like the regular solutions, 
the black hole solutions are unstable ||. 

Like EYM theory, Einstein- Yang-Mills-dilaton (EYMD) possesses static, spherically 
symmetric non-abelian regular and black hole solutions ||. Here the dilaton coupling 
constant 7 represents a parameter. In the limit 7 — > the dilaton decouples and EYM 
theory is obtained, for 7 = 1 contact with the low energy effective action of string 
theory is made, and in the limit 7 — > 00 gravity decouples and Yang-Mills-dilaton 
(YMD) theory is obtained. 

Recently we have shown, that YMD theory possesses also static axially symmetric 
regular solutions 0. These solutions are labelled by the winding number n > 1 and 
the node number k of the gauge fields. For n — 1 one obtains the static spherically 
symmetric solutions of YMD theory ||. The axially symmetric solutions have a torus- 
like shape. Choosing the z-axis as the symmetry axis, the energy density has a strong 
peak along the p-axis and decreases monotonically along the z-axis. 

The existence of the regular static axially symmetric YMD solutions strongly sug- 
gests the existence of the corresponding EYMD solutions for finite values of the dilaton 
coupling constant. In this letter we present strong (numerical) evidence, that such 
regular gravitating axially symmetric solutions indeed exist in EYMD theory and also 
in EYM theory. 

2 Axially symmetric ansatz 

We consider the SU(2) Einstein- Yang-Mills-dilaton action 

s = /(i4 + L ")^ (1 » 

with 

L M = ~d^d^ - e 2K ^Tr{F^F^) . (2) 
To obtain static axially symmetric solutions we employ isotropic coordinates and 
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adopt the metric 



fdt 2 + j (dp 2 + dz 2 ) + i 



./ 



p 2 d<p 2 , 



(3) 



with /, m and I being only functions of p and z. The corresponding ansatz for the 
purely magnetic gauge field (A = 0) is |7|, || 



n„ ,3 
1 5 



.4 
1 



(4) 

(5) 
(6) 



with the Pauli matrices r and r" = r • (cosn0, sinn0, 0), = r ■ (— sinn0, cosncf), 0). 
The four functions and the dilaton function $ depend only on p and 2;. 
With this ansatz the energy density e = — Tq = — Lu becomes 
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gw^wl + (d z wl + gw\w\ 



(7) 



The system possesses a residual abelian gauge invariance pL 0, H. With respect to 
the transformation 

U = e ir{p ' z K (8) 

the functions (pw\, pw\—njg) transform like a scalar doublet, and the functions {w\ , w^) 
transform like a 2-dimensional gauge field. We fix the gauge by choosing the gauge 
condition g, g g 



d p w{ + d z w% = 



(9) 



To make contact with the spherically symmetric case n — 1, we introduce the 
coordinates r and 8 (p = r sin9, z = r cos 9) and the gauge field functions Fi(r, 9) 0, U 



Wo 



— (1 -Fi)cosfl 
gr 

Tl 

(1 -F 3 )cos# 

gr 



Wo 



— (1-F 2 )sin0 , 
gr 

Tl 

— (1 -F 4 )sin# . 
gr 



(10) 



The spherically symmetric ansatz of ref. |4j] is recovered for Fi(r, 
F 3 (r, 9) = F 4 (r, 9) = w(r), $(r, 9) = <j>{r) and n — 1. 



F 2 {r,9) 



3 



The above ansatz and gauge choice yield a set of coupled partial differential equa- 
tions for the metric and the matter field functions. To obtain regular asymptotically 
flat solutions with finite energy density we impose at the origin (r = 0) the boundary 
conditions 

d r f = d r m = d r l = <9 r $ = , F 1 = F 2 = F 3 = F 4 = 1 , (11) 
and at infinity (r = oo) 

f = m = l = l > $ = o , F 1 = F 2 = F 3 = Fi = ±1 , (12) 

further we impose for all functions that their derivatives with respect to 9 vanish on 
the p- and the z-axis 0. The boundary conditions for the gauge field functions at 
infinity imply, that the solutions are magnetically neutral. Note, that a finite value 
of the dilaton field at infinity can always be transformed to zero via $ — ► $ — ^(oo), 

The mass M of the regular axially symmetric solutions can be obtained directly from 
the total energy-momentum "tensor" t^ u of matter and gravitation, M — j r 00 d 3 r [10] 



or from M = — J (2T ° — T^j ^—gdrdOdcf). Both expressions give the same values for 
the mass of the solutions. 

We now remove the dependence on the coupling constants k and g from the dif- 
ferential equations by changing to the dimensionless coordinate x = (g/v4irG)r, the 
dimensionless dilaton function ip = VAttGQ and the dimensionless coupling constant 
7 = k/V^ttG (7 = 1 corresponds to string theory). The dimensionless mass is then 
given by /X = (g/y/4nG)GM. 

In the spherically symmetric case the following relations between the metric and 



the dilaton field hold [11 



(p(x) = -7ln(-0 tt ) , (13) 

H = j-x 2 d x f\ oa = -x 2 d x yD\ 00 = — , (14) 
277 

where D is the dilaton charge. These relations also hold for the regular axially sym- 
metric solutions considered here. Their derivation is based on the equation of motion 



of the dilaton field and will be given elsewhere [12 . 



3 Solutions 

Subject to the above boundary conditions, we solve the equations numerically. To map 
spatial infinity to the finite value x = 1, we employ the radial coordinate x = The 
numerical calculations are based on the Newton-Raphson method. The equations are 
discretized on a non-equidistant grid in x and 9. Typical grids used have sizes 150 x 30, 
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covering the integration region < x < 1 and < 6 < ir/2. The numerical error for 
the functions is estimated to be on the order of 10~ 3 and 1CT 2 for k < 4 and k — 4, 
respectively. 

In Tables 1 and 2 we show the dimensionless mass of a subset of the regular axially 
symmetric solutions obtained so far. The energy density e = — Tq of the solutions has 
a strong peak along the p-axis, and it decreases monotonically along the z-axis. Thus 
equal density contours reveal a torus-like shape of the solutions. As a typical example 
we show the energy density e of the solution with n = 3, k = 3 and dilaton coupling 
constant 7 = 1 in Fig. |l|. 

With n and 7 fixed and increasing k, the location of the peak of the energy density 
moves inward and the peak increases in height, whereas with fixed k and 7 and increas- 
ing n the peak of the energy density moves outward and decreases in height. This is 
demonstrated in Table 2. 

The gauge field functions Fi and the dilaton function ip of the regular axially sym- 
metric EYMD and EYM solutions look similar to those of the corresponding YMD 
solutions ||. Like the dilaton function <p, the metric functions do not exhibit a strong 
angular dependence. These functions will be exhibited elsewhere [f[2| . 

For fixed n and 7, with increasing k the sequence of axially symmetric solutions 
tends to a limiting solution. The gauge field functions (Fi)k tend to the limiting func- 
tion Foo = 0. (Because of the boundary conditions at the origin and at infinity, they 
approach the limiting function nonuniformly. ) The dilaton functions tpi- tend to the 
limiting function ip^, which represents the dilaton function of the "extremal" EMD 
solution |13| with n units of magnetic charge and the same value of 7. This is demon- 



strated in Fig. Q for n = 3 and 7 = 1. We observe, that (pk deviates from ip^ only in 
an inner region, which decreases exponentially with k. 

For finite values of 7 and fixed n, the sequences of solutions thus approach as limiting 



solutions the "extremal" EMD solutions |13| with magnetic charge n and the same 7 



This generalizes the corresponding observation for the spherically symmetric solutions 



with n — 1 H, [TTJ] . For 7 = the limiting solutions represent the exterior of extremal 
Reissner-Nordstr0m black holes with charge n. For the YMD solutions the limiting 
solutions are magnetic monopoles with n units of charge ||. 



The limiting values for the mass, [l = nj ^1 + 7 2 [O, represent upper bounds 



for the sequences, as observed from Tables 1 and 2. The larger n, the slower is the 
convergence to the limiting solution. Further details on the convergence properties of 



the sequences of solutions will be given elsewhere [12 



4 Conclusions 

In addition to the known static spherically symmetric solutions, both EYMD and EYM 
theory possess sequences of regular static axially symmetric solutions. These sequences 
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are characterized by the winding number n > 1, and the solutions within each sequence 
by the node number k. (For n = 1 the spherically symmetric solutions are recovered.) 
For fixed n and 7, with increasing k the solutions tend to the "extremal" Einstein- 
Maxwell-Dilaton solution [13[ with n units of magnetic charge and the same 7. 

The multisphalerons have a torus-like shape. Apart from that, many properties 
of the axially symmetric solutions are similar to those of their spherically symmetric 
counterparts. In particular, there is all reason to believe, that these regular static 
axially symmetric EYMD and EYM solutions are unstable. Since we can also associate 
the Chern-Simons number Nqs = n/2 |§ (for odd k [Q) with these solutions, we 
interpret them as gravitating multisphalerons. 

Having constructed the axially symmetric solutions in EYMD and EYM theory, 
it appears straightforward to construct analogous solutions is theories with a Higgs 
field. We therefore expect to find gravitating axially symmetric multimonopoles (for 
the case of a Higgs triplet) or gravitating axially symmetric multisphalerons (for the 
case of a Higgs doublet). Similarly there should be gravitating axially symmetric mul- 
tiskyrmions. Work along these lines is in progress. 

We consider the above set of solutions to be the simplest type of gravitating non- 
spherical regular solutions of EYMD and EYM theory. We conjecture, that there 
are gravitating regular solutions with much more complex shapes and only discrete 
symmetries left. This conjecture is based on the observation, that for some types 
of solitons in flat space the symmetry structure of the (energetically lowest) solutions 
becomes increasingly complex with increasing winding number or charge n. For instance 
for skyrmions, the lowest n = 1 solution is spherically symmetric, the lowest n = 2 
solution has axial symmetry, and the lowest n > 3 solutions respect only discrete 



crystal-like symmetries |T5 



But EYMD and EYM theory also possess black hole solutions. The non-abelian 
spherically symmetric black hole solutions may be regarded as black holes inside sphale- 
rons [nj. We conjecture, that also the gravitating axially symmetric solutions can 
accommodate black holes in their interior. And this conjecture naturally extends to 
gravitating solutions with more complex shapes and less symmetry. The existence of 
such black hole solutions without rotational symmetry inside multimonopoles has also 
been conjectured from a stability argument [IB]. 
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EYM 


EYMD 


YMD 


fc/7 





0.5 


1.0 


2.0 


oo 


1 


1.870 


1.659 


1.297 


0.811 


1.800 


2 


2.524 


2.250 


1.770 


1.114 


2.482 


3 


2.805 


2.505 


1.976 


1.247 


2.785 


4 


2.922 


2.611 


2.063 


1.304 


2.913 




3 


2.683 


2.121 


1.342 


3 



Table 1 

The dimensionless mass fi of the EYMD solutions with winding number n = 3 and 
up to 4 nodes for several values of the dilaton coupling constant 7. For comparison, 
the last row gives the mass of the limiting solutions, the first column gives the mass 
of the EYM solutions, and the last column the scaled mass of the corresponding YMD 
solutions ||. 





fi 


k/n 


1 


2 


3 


4 


1 


0.577 


0.961 


1.297 


1.607 


2 


0.685 


1.262 


1.770 


2.239 


3 


0.703 


1.365 


1.976 


2.549 


4 


0.707 


1.399 


2.063 


2.698 




0.707 


1.414 


2.121 


2.828 






^max 


(Pmax ) 




k/n 


1 


2 


3 


4 


1 


1.075 (0.) 


0.177 (0.90) 


0.098 (1.59) 


0.072 (2.37) 


2 


11.63 (0.) 


0.910 (0.30) 


0.380 (0.66) 


0.235 (1.10) 


3 


79.70 (0.) 


3.443 (0.09) 


1.124 (0.28) 


0.601 (0.51) 


4 


498.2 (0.) 


12.01 (0.03) 


3.064 (0.11) 


1.435 (0.24) 



Table 2 

The dimensionless mass fi, the maximum of the energy density e max and the location 
Pmax of the maximum of the energy density of the EYMD solutions of the sequences 
n = 1 — 4 with node numbers k = 1 — 4 and 7 = 1. For comparison, the mass of the 
limiting solutions is also shown. Note, that fi/n decreases with n for fixed finite k. 
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6 (p , z ) n = 3 k = 3 7=1 




Figure 1: The energy density e = — Tq is shown as a function of the dimensionless 
coordinates p and z for the EYMD solution with winding number n — 3, node number 
fc = 3 and 7 = 1. 
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Dilaton Function cp k (x,9) 
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Figure 2: The dilaton functions ifik(x, 9) for the EYMD solutions with winding number 
n — 3, node numbers k — 1 — 4 and 7=1 are shown as a function of the dimensionless 
coordinate x. The dashed, the dotted and the dash-dotted lines represent the angles 
9 = 0, 9 = 7r/4 and 9 = ir/2, respectively. Also shown is the limiting function y?oo(:r) 
(solid line). 
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